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Abstract — A combined base station association and power con- 
trol problem is studied for the uplink of multichannel multicell 
cellular networks, in which each channel is used by exactly 
one cell (i.e., base station). A distributed association and power 
update algorithm is proposed and shown to converge to a Nash 
equilibrium of a noncooperative game. We consider network 
models with discrete mobiles (yielding an atomic congestion 
game), as well as a continuum of mobiles (yielding a population 
game). We find that the equilibria need not be Pareto efficient, 
nor need they be system optimal. To address the lack of system 
optimality, we propose pricing mechanisms. It is shown that these 
mechanisms can be implemented in a distributed fashion. 



I. Introduction 

Wireless communication systems have experienced tremen- 
dous growth over the last decade, and this growth continues 
unabated worldwide. Efficient management of resources is 
essential for the success of wireless cellular systems. In a 
mobile cellular system, mobiles adapt to time varying radio 
channels by adjusting base station (BS) associations and 
by controlling transmitter powers. Doing so, they not only 
maintain their quality of service (QoS) but also enhance their 
transmitters' battery lives. In addition, such controls reduce the 
network interference, thus maximizing spatial spectrum reuse. 
Distributed control is of special interest, since the alternative 
of centrally orchestrated control involves added infrastructure, 
the need for distribution of measurements, and hence system 
complexity. 

Distributed control algorithms for single channel multicell 
networks have been extensively studied (Foschini & Mil- 
janic Yates J2], Hanly O). The monograph by Chiang 
et al. (4) and references therein provide an excellent survey 
of the area. Noncooperative games have been a natural tool 
for analysis and design of distributed power control algo- 
rithms. Scutari et al. (5) and Heikkinen (6) model distributed 
power control problems as potential games, while Altman & 
Altman (7) show that many of the cellular power control 
algorithms can be modeled as submodular games. In contrast, 
uplink resource allocation for multichannel multicell networks 
poses several challenges as observed in Yates (2 and Jiang 
et al. 0. 
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We address the resource allocation problem in the uplink 
of a multichannel multicell network with a single traffic class. 
Such a problem arises when a CDMA operator chooses to 
lease and utilize multiple frequency bands in order to reduce 
in-network interference, or multiple operators who lease differ- 
ent bands decide to cooperate. The newer mobile devices are 
typically radio agile, and thus have the option to choose from 
one of these distinct bands. We address a simplified version 
of this multichannel multicell problem where all BSs operate 
on different frequency bands. 

A preview of our results is as follows. We propose a 
distributed algorithm for the combined base station association 
and power control problem, and subsequently model the 
problem as a player-specific congestion game. The equilibrium 
states of such algorithms, which are Nash equilibria of the 
corresponding games, may be far from system optimum. We 
resort to pricing mechanisms to induce mobiles to behave in 
a way that optimizes system cost. We also show that such a 
mechanism can be employed in a distributed fashion. Towards 
this end, we model the network as having a continuum 
of (nonatomic) mobiles, each offering infinitesimal load, which 
leads to a population game formulation. We then provide a 
marginal pricing mechanism that motivates a pricing strategy 
for the discrete mobiles case. Note that, unlike the case 
of transportation networks, mobiles are not really priced in 
cellular networks. The pricing is simply a part of the decision 
making routine built into each mobile in order bring about a 
distributed control mechanism that drives the system toward 
optimality. 

The paper is organized as follows. In Section ITU we briefly 
discuss concepts of finite noncooperative games and popula- 
tion games. We study a network model with discrete mobiles 
in Section [HI] We propose a combined association and power 
control algorithm, model it as a noncooperative game, and 
analyze its performance. We extend this analysis to a network 
with a continuum of mobiles in Section [IV] To address 
the inefficiency of the proposed algorithms, we design toll 
mechanisms in Section [V] Finally, we conclude the paper with 
some remarks in Section [VT] 

Optimal power allocation and BS association in uplinks 
of multichannel multicell cellular networks have not been 
explored before. Ours is an attempt at a detailed coverage 
on what is possible in general, with more specific results in 
some special cases. 
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II. Game Preliminaries 
A. Finite Noncooperative Games 

A noncooperative strategic form game ( Ai , (Ai , i £ 
Ai),(ci,i £ -M)) consists of a set of players Ai = 
{1, ...,M}. Each player i is accompanied by an action 
set Ai and a cost function : X-fLiAi — > R. In this 
work, we assume all action sets to be finite. An action 
profile a = (a,i,i = 1,...,M) prescribes an action <Zj 
for every player i £ Ai. For a = (a.j,2 = 1,...,M), 
denote a , := (ai, . . . , aj_i, ai+i, . . . , om) and (6i,a_i) := 
(oi, . . . , ai_i, 6j, a»+i, . . . , aA/). 

Definition 2.1: Nash Equilibrium (NE): For an action pro- 
file a, a mobile i's best response, 2?,: (a) C Ai, is defined 
as Sj(a) := argmin;,. a(bi, a_i). a is said to be a Nash 
Equilibrium for the game if ai £ Si (a) for all i g AL 

Definition 2.2: Potential Game: A game (Ai, (Ai,i £ 
Ai), (ci,i £ Ai)) is said to be an ordinal potential game if 
there exists a function V : ~xf£iAi —> R, known as an ordinal 
potential function, that satisfies 

Cj(6i,a_i) < c»(a) ^ V(6*,a_i) < V(a) 

for all i £ Ai, h £ Ai, a £ x^A*. 

Clearly, all minimizers of an ordinal potential function V 
are Nash equilibria of the game. Thus all ordinal potential 
games (Ai, (Ai,i £ Ai), (cj, i £ Ai)) admit at least one Nash 
equilibrium. On account of their finiteness they also have the 
finite improvement path (FIP) property, i.e., do not contain 
improvement cycles (Monderer & Shapley J5] Lemma 2.3]). 
Thus, in a finite ordinal potential game when players update 
as per the better response strategy, round-robin or random 
update processes converge to a Nash equilibrium in a finite 
number of steps. With the same strategies, an asynchronous 
update process also converges (Neel ifTol Chapter 5]). 

Remark 2.1 : The strategic form games that have the FIP 
property also admit the finite best-response path (FBRP) prop- 
erty, i.e., do not contain best response cycles (Milchtaich ifTTl 
Section 5])Q Thus, if players update as per the best response 
strategy, then also the above update processes converge to 
a Nash equilibrium in a finite number of steps. The reverse 
implication is not true in general - FBRP need not imply FIP. 

Definition 2.3: Congestion Game: A game (Ai, (Ai,i £ 
Ai), (ci,i £ Ai)) is said to be a player-specific weighted 
singleton congestion game if 

1) there exists a set AS such that Ai = Af for all i £ Ai, 
and 

2) there exist constants £ Ai) and nonincreasing 
functions /y , i £ Ai, j £ Af such that 

Ci(a) = / itti ( ^2 Pi) for a11 a G x tUA,i £ Ai. 

ieM: 
ai —ai 

In the above definition, we interpret Af as a set of facilities and 
Bi as the load offered by player i. Then, Y] ieM . Bi denotes the 

ai=at 

total load on facility a,, under an action profile a. The game is 

1 A best response cycle is a finite best response path a 1 , ■ ■ ■ , a k such that 
a 1 = a k , and for some j S {1, • • ■ , k — 1}, the deviating player in iteration 
j strictly benefits ll2l . 



a singleton congestion game because each action picks exactly 
one facility. It is weighted because players offer different loads, 
and it is player-specific because the cost functions c,(-) are 
player-specific. 

Rosenthal lfT3l has defined congestion games with un- 
weighted players and player-independent cost functions, but 
more general action sets. The above generalization is due to 
Milchtaich IfTTl who also showed that singleton congestion 
games with weighted players but player independent costs 
admit FIP. Gairing et al. lfT4l studied these games in the 
special case of linear cost functions. Georgiou et al. lfT31 
introduced a routing game where the players only have incom- 
plete information about edge latencies, and showed that this 
game can be transformed into a (complete information) player 
specific singleton congestion game with linear cost functions. 
Mavronicolas et al. Ifl6l considered a subclass where each 
player-specific cost function is composed (by means of an 
abelian group operation) of a player-specific constant and a 
facility-specific nondecreasing function. 

B. Population Games 

A population game (Sandholm IfTTl ) (Ai,(Ai,l £ 
£),(cij,l £ C,j £ Ai)) consists of C = {1,...,L} classes 
of nonatomic populations of players. Ai = U; £ £7\4;, and 
Mi := \Aii \ denotes the total mass of the class / population. 
By a nonatomic population, we mean that the mass of each 
member of the population is infinitesimal. Players of class I 
are associated with an action set Ai- Actions of these (class 
I) players lead to an action distribution m' = (mij,j £ Ai), 
where 2~2jeAi mi i = ^ me pl avers within a class are 
alike. Thus the action distributions completely specify the 
play; we can characterize the states and dynamics of play 
solely in terms of action distributions. Let m = (m , I £ C) 
denote the action distribution profile across the entire pop- 
ulation, and let Ai* denote the set of all such profiles. A 
population I is also accompanied by continuous cost density 
functions cy : Ai* — > R. 

Definition 2.4: Nash Equilibrium (NE): An action distribu- 
tion profile m is a pure strategy Nash equilibrium for the 
game (Ai, (A,,i £ Ai), (cy,Z £ C,j £ (Ai)) if and only if 
for all I £ C and j £ Ai, a positive mass my > implies 
cy(m) < Qfc(m) for all k £ Ai. 

Remark 2.2: At a Nash equilibrium m, for a class I, if j 
and k are any two facilities in Ai such that my > 0, mik > 0, 
then cy(m) = Q fc (m). 

Definition 2.5: Potential Game: A game (Ai,(Ai,l £ 
C), (cij, I £ C,j £ Ai)) is said to be a potential game if there 
exists a C 1 function V : Ai* — > R, known as a potential 
function, that satisfies 



for all I £ C,j £ Ai,m £ Ai*. 

It is well known that Nash equilibria are the profiles which 
satisfy the Kuhn-Tucker first order conditions for a minimizer 
of the potential function (Sandholm IfTTl Proposition 3.1]). Any 
dynamics with positive correlation and noncomplacency (in 
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particular the best response dynamics) approaches a Nash 
equilibrium ifTTl . 

We are interested in nonatomic congestion games (Sand- 
holm Q7)), in which Ai = Af,Vl e £, for a given set TV. 
As before, we interpret Af as a set of facilities. Moreover, 
each class / has an associated offered load density ji > 0. 
An action distribution profile m leads to a congestion profile 
(irij,j e Af), where rrij = J^;e£ m y7(- The cost density 
functions cy depend on m only through rrij, and are increasing 
in rrij. 

C. Pricing 

Levying of tolls is a conventional way to enforce system 
optimality in nonatomic networks. Beckman |[T8l and Dafer- 
mos & Sparrow |fT9l studied optimal tolls in transportation 
networks with a single class of users. Later Dafermos 11201 
and Smith |2T1 extended the analysis to multiclass networks. 
Roughgarden & Tardos |22| applied these ideas in computer 
networks and analyzed tolls for optimal routing. 

In the atomic (discrete) setting, Caragiannis et al. fl23l 
proposed tolls for two-terminal parallel-edge networks with 
unweighted users and linear latency functions. Fotakis & 
Spirakis 11241 extended these to generic two-terminal net- 
works with unweighted users and arbitrary increasing latency 
functions. We propose an alternative toll mechanism, and 
demonstrate that the proposed tolls can be computed in a 
distributed fashion. 

III. Discrete Mobiles 

A. System Model 

We now describe the model adopted in this work. We 
consider the uplink of a cellular network consisting of several 
BSs and mobiles. Each BS operates in a distinct frequency 
band. Let Af = {1, . . . , N} and M = {1, . . . , M} denote the 
set of BSs and the set of mobiles, respectively. 

A mobile must be associated with one BS at any time, and 
is free to choose the BS with which it associates. Let hij 
denote the power gain from mobile i to BS j. Let the receiver 
noise at all BSs have the average power a 2 . Let Pi denote the 
power transmitted by mobile i, and let Oj be the BS to which 
it is associated. Under an association profile a = (di,i = 
1, . . . , M), let A4j{&) be the set of mobiles associated with 
BS j. Under an association profile a = (<Zj,« = 1,...,M) 
and a power vector p = (pi,i = 1, ...,M), the signal to 
interference plus noise ratio (SINR) of mobile i at BS a, is 

hiaiPi 

J2l£M az {a)\{i} hlo-tPl + C 2 

Mobile i has a target SINR requirement ji. 

B. The Proposed Algorithm 

Yates (2j and Hanly Q proposed an algorithm for dis- 
tributed association and power control in single channel cel- 
lular networks. Convergence results for the algorithm are 
based on the concept of a standard interference function. The 
technique is based on a mobile reassociating itself with a 



BS with which it needs to use the least power; this fails to 
work in the case of a multichannel network and analogous 
convergence results for this algorithm may not hold (see 
Yates J2] Section VI]). Even in instances where the algorithm 
converges, it may get stuck at a power allocation that is not 
Pareto efficient. 

We propose an alternative distributed algorithm for com- 
bined BS association and power control in multichannel mul- 
ticell cellular networks. We also show its convergence. We 
make use of the following simple fact (see, for example, 
Kumar et al. |25l Chapter 5]). Consider the subproblem of 
power control with a fixed association a. Define /3j = jfz^, a 
measure of the "load" offered by mobile I. 

Proposition 3.1: For a fixed association a, 

(i) The power control subproblem of BS j is feasible iff 

E(£JMi(a) A < 1- 

(ii) If the power control subproblem of BS j is feasible, there 
exists a unique Pareto efficient power vector p(a) given by 

pM) = -j—-, — ^ o- 

Throughout we assume that there exists at least one feasible 
association and power vector. Proposition 13.11 motivates the 
following algorithm. 

Multichannel Association and Power Control (MAPC): 

Mobiles switch associations in a round-robin fashion by taking 
into account the optimal power consumptions (given by Propo- 
sition |3.1f ii)) at the BSs with which these associate. This is 
done as follows. As the load at a BS changes, it immediately 
broadcasts the new load, and associated mobiles update their 
powers to the optimal required powers as per the new loads. 
To be more precise, define 

where [x} + = max(x, 0). For t = 0, 1, 2, . . . , mobile i where 
i = l + (t mod M) updates its association and power at t+l 
if a,i(t) £ argmin J - gAr c i ((i,a(t)_ i )). In this case, 

a,(t + 1) G argminci((i,a(t)_i)), (2a) 

and with a(t + 1) = (a^t + l),a(t)_i), 
Pl(t + 1) =q(a(t + l)), 

V/ G M 0i(t) (a(i)) U M adt+1) (a(t + 1)). (2b) 

Remark 3.1: A mobile i should not choose a BS if the 
device renders the corresponding power control subproblem 
infeasible. The situation is characterized by YlieM (a) @ l — 
and <[TJ appropriately yields infinite cost for the mobile. Even 
if the algorithm starts with an infeasible association, selfish 
moves of players eventually lead to a feasible one, and updates 
remain feasible thereafter. 

Note that while only one mobile updates its association 
at a time, all mobiles that perceive a change in load at 
their BSs update their powers to optimal values based on the 
new loads. Simultaneous association updates are not allowed. 
In a framework with no synchronizing agent and with an 
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arbitrarily fine time-scale, it is unlikely that two mobiles 
update simultaneously. If two or more BSs result in the same 
steady state power, one is chosen at random by the mobile. 

This algorithm is also distributed in nature as the one 
proposed in Yates (2j. BS j broadcasts its total congestion 
SzeA-ij(a) ft- I R addition, each mobile i is told its scaled gains 
^ by each BS j G TV. 

C. A Congestion Game Formulation 

To show the convergence properties of the proposed algo- 
rithm, we model the system as a strategic form game. Let the 
mobiles be the players and the action set for each player be 
the possible associations, i.e, A% = TV for all i G M. Define 
the cost functions of the players to be c;(a) for all i G M. It 
can be seen that above is a player-specific singleton weighted 
congestion game, and belongs to the subclass of congestion 
games with multiplicative player-specific constants described 
in lfl6l . In the following we refer to it as the strategic form 
game (M,Af, (c l: i G M)). 

Before analyzing the general game, we consider the follow- 
ing special cases. 

1) Single Class Traffic: This is the case where all the mo- 
biles have a common target SINR requirement 7; /3 := j^—- 
In this case, 



Ci(a) 



ft 



h iai [1 - \M at (s)\ft] + 

and we have a player specific unweighted singleton congestion 
game. 

2) Collocated Users: In this case, all mobiles are situated 
close together in a group. Thus hij = hj for all i G M, j G TV, 
and 



Ci(a) = 



ft, 



V [l-EiGA-l ai (a)A] + ' 



This yields a player independent weighted singleton conges- 
tion game. 

3) Collocated BS: Here all BSs are assumed to be situated 
close together. Thus hij = hi for all i G M , j G TV, and 

Ci(a) 



ft 



hi [l-EteM^ft^' 

Now, we get a player specific weighted singleton congestion 
game. 

The following result ensures that MAPC converges in each 
of these special cases. 

Proposition 3.2: The finite strategic form game 
(.MjTV, (ci,i G .M)) is an ordinal potential game and 
thus admits the FIP in each of the following cases. 

1) = f3 for all i G M, 

2) hij = hj for all i G M , j G TV, 

3) h^ = hi for all i G M,j GTV. 

Proof: In each case, we show that the game 
(M,Af, (ci,i € M.)) is better response equivalent (Neel IflOl 
Chapter 5]) to an ordinal potential game (by demonstrating an 
ordinal potential function for the latter). This implies that, in 
each case, (Ai,Af, (ci,i £ -M)) itself is an ordinal potential 



game. It is also finite which implies that the FIP property 
holds. 

1) The strategic form game (Ai , TV, (c, , i G M.)) is better 
response equivalent to (Ai,Af, G -M)). Also note that 



1 



[l-|M ai (a)|/3] 



Ci(a) a 2 ft 
The function V 1 : N M R given by 

(\M k (B)\ 

-^n^n n ^ 

' leM keM \ t=i 



satisfies 

Vi&a-O-VxCa) 



1 



1 



Cj(j,a_i) Q(a) 



n ^ 

ievw\{i} 



/|M fc (a)\{i}| 

xll n i 1 -^" 

for all j G M,j G TV, a G TV^. Thus the game 
(A^,TV, (— — , i E is an ordinal potential game with a 

potential function V^i0 

2) The strategic form game (A4,Af,(ci,i G is better 

response equivalent to (A^,TV, (— — ,i G A4)). Also note that 



ft 



;e» ai ( a ) ft] + 



Ci(a) 

For the function V" 2 : N M 



Va(a) 




F 2 (j,a_ i )-F 2 (a) = -2 



given by 

E a 



(i - fti) 
ft 



,Ci(j,a-j) Ci(a), 

for all i G TW,j G TV, a G TV^. Thus Va is a potential 
function for the game (M,J\f,(-§,i G M)), and so the 
latter is an ordinal potential game. 

3) The strategic form game (A4,Af,(ci,i G M)) is better 
response equivalent to (TW,TV, (— 7^-, « G .M)). Also note 



that 



ft 



/ijC^a) 
The function F 3 : TV 7 ^ ->• I 



[1-E^(a)ft] + 

a 2 

M defined as 

ft[i-E ieMo .ft] + 



satisfies 

F 3 (i J a_ i )-F 3 (a) = -2^ 



ft 



ft 



hiCi(j,a^i) h i Ci(a) / 

for all z G .M,j G TV, a G J\f M . Thus the game 
(M. , TV , (— , i G .M)) is an ordinal potential game with 
V3 as a potential function. ■ 

2 This potential function is similar to those proposed in [14] for linear cost 
functions, and in 1161 for cost functions composed of player-specific constants 
and facility-specific functions. 
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Now, we focus on the general case. Gairing et al. T\4j 
show (via a counter-example with 3 players) that player- 
specific weighted singleton congestion games with linear cost 
functions are not necessarily ordinal potential games, and so, 
need not possess FIR This negative result applies to our game 
also, and convergence proofs based on potential functions 
cannot be used. However, it follows from IfTTI that the strategic 
form game (Ai,AT,(ci,i £ -M)) admits (i) FIP property if 
\Af\ = 2, (ii) FBRP property if \M\ = 2. 

Georgiou et al. |[T5l establish that player-specific weighted 
singleton congestion games with 3 players and linear cost 
functions possess FBRP property!^ Specifically, they show in 
an exhaustive manner that such games do not possess any best 
response cycles^ Their result and proof technique extend to 
the game (M,Af, (c%,i £ M.)) even though the cost functions 
Ci are not linear. Thus, the game (A4,Af, (ci,i £ -M.)) can be 
shown to possess FBRP if \M\ = 3. 

In case of more than 3 players, convergence of the best 
response dynamics in weighted singleton congestion games 
with linear cost functions is an open problem lfT31 . l26l . 
Georgiou et al. [fT31 conjecture that such games always admit 
at least one NE. Though functions Cj are not linear, the game 
(A4,Af,(ci,i £ M.)) is best response equivalent to another 
game in which costs are composed of multiplicative player- 
specific constants and affine nondecreasing functions. Also, 
simulations run on numerous instances of the game suggest 
that players' updates as per the best response strategy always 
converge in a finite number of steps. We therefore conjecture 
that 

Conjecture 3.1: The finite strategic form game 
(M, Af, (ci, i £ M.)) admits FBRP and thus possesses 
at least one pure strategy Nash equilibrium. 

The FBRP property ensures that MAPC converges in a fi- 
nite number of steps (see Section lH-Al i. Consider the following 
variants of MAPC. 

1) At each t, one mobile is randomly chosen to update its 
association. All mobiles have strictly positive probabili- 
ties of being chosen. 

2) At each t, each mobile i updates its association with 
probability e.; £ (0, 1). There is thus a strictly positive 
probability that any subset of mobiles may update their 
associations simultaneously. As before, all mobiles up- 
date their powers based on the new loads. This algorithm 
does not require any coordination among mobiles (to 
ensure one by one updates), and is thus fully distributed. 

The FBRP property of the game (M,J\f,(ci,i £ M)) im- 
plies that these two algorithms also converge to a NE with 
probability 1 (see Neel |10| Chapter 5]). 

D. System Optimality 

A system optimal power allocation should bring about the 
lowest interference environment. This motivates the following 

3 On the other hand, Mavronicolas et al. 1161 demonstrate a best response 
cycle in a game with 3 players and costs composed of additive player-specific 
constants and facility-specific nondecreasing functions. 

4 There does not seem to be any reason why this technique cannot be 
extended to more than 3 players; though the number of possibilities in the 
exhaustive search may become enormous. 



definition of system optimality. 

Definition 3.1: For an association profile a, define a system 
performance measure C(a) = X]fc-i c i( a ) w i tn c i( a ) defined 
in ((T). We define an association profile a° to be system 
optimal if it minimizes C(a) over all possible associations 

Let us now recall the following notion of Pareto effi- 
ciency j25l Chapter 5]. 

Definition 3.2: An association profile a is said to be Pareto 
dominated by another association profile a' if Cj(a') < Cj(a) 
for all i £ A4 with Cj(a') < Cj(a) for some i. An association 
profile a is said to be Pareto efficient if it is not Pareto 
dominated by any other association profile in x f£ lw 4,;. 

Clearly any association profile that is system optimal is 
also Pareto efficient. Thus, if there is a unique Pareto efficient 
association profile, it is also the unique system optimal one. 
However, unlike the case of single channel networks, joint 
association and power control problems in multichannel net- 
works do not in general admit a unique Pareto efficient power 
allocation. In particular, when > |7V|, there cannot be 
unique Pareto efficient power allocation]! To see this, define 
<di for any mobile i as the set of best match BSs as follows 

c>, := argmin— — 

The system optimal association profile a° is clearly Pareto 
efficient. Next, two cases are possible. 

1) For all i, a° £ 9,. Since \M\ > \Af\, there exist two 
mobiles i and I such that a° = af. 

2) There exists a mobile i such that a° ^ 0,. 

Consider a mobile i as in Case 1, or as in Case 2. Let a' 
be another profile which is system optimal subject to i being 
associated with any of its best match BSs and no other mobile 
being associated with that BS. It can be easily checked that 
a' is also Pareto efficient. 

As the following example illustrates, MAPC may settle at 
a Pareto inefficient association profile, and hence may not be 
system optimal. 

Example 3.1: Consider a network with two BSs, two mo- 
biles, and a common SINR requirement 7. The two BSs 
operate in disjoint bands. Assume 

hi2 < /in < 7- r 

(I-7) 

and /121 < "22 < 7 r- 

(I-7) 

The unique Pareto efficient association is (ai = 1,<22 = 2) 
with power allocation (7^777, 7^777)- However, if we start with 
initial association (<Xi = 2,a 2 = 1), MAPC will not move 
forward, because a unilateral switch requires larger power to 
meet the target SINR. Neither mobile will switch to the BS 

2 2 

with which it has a better channel. Hence, (7^7-7, 7^777) i s a 
steady state power vector at which the algorithm settles; it is 
Pareto inefficient. 

5 Though two Pareto efficient power allocations can be identical up to a 
permutation, e.g., if two mobiles are indifferent with respect to their SINR 
requirements and channel gains to all the BSs. 
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In the following we consider special cases, and investigate 
whether the proposed algorithm leads to a system optimal 
association profile. 

1) Collocated Mobiles and Single Class Traffic: Even in 
this special case, MAPC may settle at a Pareto inefficient NE 
as shown in the following example. 

Example 3.2: Consider a 2-cell network with 4 collocated 
mobiles and = f3,i = 1,2,3,4. Assume that hi and hi 
satisfy 



Mi 
Mi 



3/3) 
2/3) 



> 



Mi 
Mi 



2/3), 
/3). 



The following facts are easily verified. Both the inequalities 
can be met simultaneously. The association {a\ = a 2 = CI3 
1 . 04 = 2) is a NE from which the algorithm does not move. 
This association is Pareto dominated by (03 = 04 = 1, 01 = 
02 = 2) which is another NE. Thus MAPC may settle at a 
Pareto inefficient NE. 

Consider now a variant of MAPC in which mobile i = l + (t 
mod M) updates its association at i + 1 if] 



a,i(t) g argmin a(j,a(t). 



-i), E # 

leMjUrtt)-*) 



In this case, 



Oi(t + 1) G argmin | Cj(j, a(t)_j), 



E 



In words, a mobiles selects a least loaded BS (counting the 
tagged mobile also) among the ones which require transmis- 
sion with the least power. We name this variant MAPC*. 

Proposition 3.3: If the strategic form game (A4,M, (ci,i G 
Mi)) contains no best response cycles, then MAPC* converges 
in a finite number of steps. 

Proof: Suppose a° is the initial association profile. Sup- 
pose a 1 , • ■ ■ , a k are successive association profiles generated 
by MAPC*, with a 1 possibly the initial association profile 
a°. For an association profile a, let m(a) be the congestion 
vector with its elements arranged in a decreasing order. The 
following two cases are possible. 

1) For some j G {1, ■■■ ,k — 1}, the deviating mobile 
strictly benefits. Then, a k 7^ a 1 because the game 
(A4,M, (ci,i G Ml)) does not contain any best response 
cycle. 

2) The deviating mobiles do not strictly benefit for any j G 
{1, ■ ■ • , k — 1}. Then, it must be that m(a 1 ), ■ • • , m(a k ) 
is a lexicographically decreasing sequence. Again, a k 7^ 
a 1 . 

We thus conclude that iterates generated by MAPC* never 
contain a cycle. Since there is only a finite number of distinct 
association profiles, MAPC* must converge in a finite number 
of steps. ■ 
The FBRP property (Conjecture [Til ensures that MAPC* 
converges in a finite number of steps. We now show that 

6 The minimization is with respect to the lexicographical ordering. 



MAPC* converges to a Pareto efficient NE in the special case 
of collocated mobiles and single class traffic. 

Proposition 3.4: For the noncooperative game 
(M.,Af,(ci,i G M)), when M = hj an( l A — P f° r 
all i G M, j G M, the steady states of MAPC* are the Pareto 
efficient NEs of the game. 

Proof: Propositions 13.21 and 13.31 imply that MAPC* con- 
verges in a finite number of steps in this special case. For 
any association profile a, let rrij (a) be the number of mobiles 
associated with BS j. Let a be a NE, and a' be another profile 
dominating a. We show that the proposed variant of MAPC 
does not settle at a. 

We first argue that congestion vectors m(a) = 
(mjv(a), • • • , mjv(a)) and m(a') cannot be identical. Indeed 
if this is the case, a' is obtained by permuting the mobiles' 
associations in a in some way. But then their payoffs undergo 
the same permutation, which makes it impossible for all of 
them to gain. 

We define g.j = j- and /(to) = rzz§^m+ ■ Then a mobile 
associated with BS j incurs a cost gjf(rrij(a)). Further, a 
being a NE, 

ft/K'(a)) < gkf{m k (a) + 1) 

for all j, k G M. In particular, 

TOj(a) < TOfc(a) + 1 if g k < g j} (3a) 
and TOj(a) < m k (a) + 1 if g k < gj. (3b) 

Next, we define 

c := .max 3j/K(a)), 

jGN-.mj (a)>0 

and Af\ := argmax gjf(mj(a)) 

j£A/':m i (a)>0 

Under a' none of the mobiles incurs a cost more than c. In 
particular, those associated with a BS j G A/i under a' must 
have cost less than c. This implies m,j(a') < m,j(a) for all 



j G Afi- Now suppose that rrij(a') 



(a) for all j G Mi, 



and mfc(a') > m k (a) for a k G M \ Mi. Then, 

g k f(m k (a')) > g k f(m k (a) + 1) > gjf(mj(a)) 

for any j G Mi. The last inequality holds because a is an NE. 
Thus we have that 

g k f(m k (a')) > c, 

and hence there are more mobiles incurring costs greater 
than or equal to c under a' than under a. This contradicts 
the hypothesis that a' Pareto dominates a. Thus there must 

be BSs j G Mi,k G M \ Mi with rrij(a') < rrij(a) and 
TOfc(a') > TOfc(a) which is same as m k (a') > m k (a) + 1. 
Again, a being an NE, 

g k f(m k (a')) > c. 

But the hypothesis that a' Pareto dominates a implies that 

g k f{m k {^)) < c. 
Thus k must belong to the set 

M 2 := {k G M \ Mi : g k f(m k (a) + 1) = c} 
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Moreover, mfe(a') = mfe(a) + 1. 

Now, we claim that there exist BSs j G Afi and fc € A/a 
such that g., < g k . Assume this claim holds. Then, 

g k f(m k (a) + 1) = c= gjf(mj(a)) 

implies that mj(a) > mfe(a) + 1. Thus, under the proposed 
algorithm, one of the mobiles associated with BS j moves to 
BS fc, i.e., the algorithm does not settle at a. 

We prove the claim via contradiction. Suppose gj > g k for 
all j G M\ , fc G A<2 ■ Obtaining a' from a may involve three 
types of load transfers. 

1) One mobile moves from a BS j G Af\ to a BS k G A/2 
such that <7j = g k . By the definition of A/2, such moves 
only permute the overall cost profile, and by themselves 
cannot lead to a'. 

2) One mobile moves from a BS j G Af\ to a BS fc G A/2 
such that gj > g k . Then, the cost reduces for rrij(a) — 1 
mobiles that are still with BS j, but increases to c for 
TOfc(a) > rrij (a) — 1 mobiles (see d3a|i). Such moves also 
cannot lead to the association profile a'. 

3) n > 1 mobiles move from a BS j G Afi to BSs 
ki, . . . ,k n G A/2 (they have to move to different BSs, 
again by the definition of A/2). Now, the cost reduces for 
rrij (a) — n mobiles, but increases to c for 

n 

V^mjfc^a) > n(jnij{a) — 1) > n%j{a) — n (see (l3bl ) 
;=i 

mobiles. Such moves also cannot lead to the association 
profile a'. 

Thus there must be BSs j G Afi and fc G A/2 such that g 3 < % 
as claimed. This completes the proof of the proposition. ■ 

However, the obtained Pareto efficient association profile 
need not be system optimal. This is demonstrated by Exam- 
ple [4j2] for the case of a continuum of mobiles. 

2) Collocated BSs: Next, we consider the case where the 
BSs are collocated, so that hij = hi for all i G M,j G Af. 
Mobiles may have different target SINR requirements. For any 
association profile a, define its support S a to be the set {j G 
Af : ai — j for some I G M}. We say that a has full support 

if 5 a = Af. 

Lemma 3.1: In the game (A4,Af, (q, i G A4)), when hij = 
hi for all i G A4,j G Af, any association profile with full 
support is Pareto efficient. 

Proof: Consider a, an association profile with full sup- 
port. Let there be another association profile a' that dominates 
a. Assume that m(a) and m(a') are the congestion vectors 
corresponding to a and a' respectively. Identify Ji , • • • , i/y € 
Af and k\, . . . , fcjv G Af such that rrij 1 (a) < • • • < m,j N (a) 
and TOfc 1 (a / ) < ■■• < nik N (a'). Then we can show by 
induction that for every I = 1, . . . , N we must have rrik, (a') < 
mj,(a). Since a' dominates a, for a mobile i with aj = ji 
and = fc r , its costs under a and a' must satisfy 

a 2 fa > o 2 Pi 

hi [1 - m jl (a)}+ hi [1 - mfc r (a')]+ ' 

But mfc r (a') > TOfc^a'), implying that 

a 2 Pi > o 2 Pi 

hi [1 - m J1 (a)}+ ~ hi [1 - m kl (a')] + ' 



Moreover, our feasibility assumption ensures that both the 
sides in the above inequality are finite. Clearly, mfc^a') < 
m,j 1 (a). Assume m.fc i _ 1 (a / ) < mj l l (a). Indeed, if m kl (a') 
were greater than m,j l (a), then in order to keep their costs at 
most as per a, all mobiles associated with the BSs ji, ■ ■ ■ ,ji as 
per a must be associated with the BSs k\, . . . , fc;-i in a'. But 
by the full support hypothesis rrij l (a) > 0, which when put 
together with the induction hypothesis yields 2~2r=i m >( a ) > 
E r =i m /!r( a ')' S° a ll °f tne mobiles associated with BSs 
ji, . . . ,ji under a cannot be associated with BSs k\, . . . , fc;-i 
under a'. Thus the induction step is proved. Next, since 

Eili m j ; ( a ) = Eili m fci( a ')' we must have ™fci( a ') = 
rrij^a) for all / = 1,...,N, But under this condition, all 
mobiles have the same cost under a and a', contradicting 
Pareto domination. ■ 
Proposition 3.5: All the Nash equilibria in the game 
(M,Af, (ci,i G M)), when hij = h t for all i G M,j G Af, 
are Pareto efficient. 

Proof: Let a* be a Nash equilibrium. The following are 
the two possible scenarios. 

1) a* does not have full support: We must have | | < \Af\. 
Indeed, if \A4\ > \Af\ and a* does not have full support, 
then there must be mobiles i and / with a* = a ; * and 
a BS j with A^ (a*) = 0. Clearly, mobile i benefits by 
moving to BS j. This contradicts the fact that a* is a NE. 
Next, i 7^ / implies a* 7^ for the same reason as 
explained above. Since all BSs have the same channel 
gain to a mobile, a* is Pareto efficient. 

2) a* has full support: Lemma IXTl implies that a* is Pareto 
efficient 

■ 

However, a NE need not be system optimal if traffic is not 
of single class or mobiles are not collocated as shown in the 
following examples. 

Example 3.3: Consider a 2-cell network with collocated 
BSs. The 4 mobiles are symmetrically located with respect 
to the 2 collocated BSs. Thus hj = h,j = 1,2. For this 
example assume Pi = 2/3, i = 1,2, and Pi = 3/3, i = 3,4. It 
can be seen that (ai = a-i = 1, 03 = 04 = 2) is a NE which 
is Pareto efficient. But its social cost is higher than that of 
[a\ = 03 = 1, a-2 = (H = 2) which is another NE. 

Example 3.4: Consider a 2-cell network with 5 mobiles. 
The 2 BSs are collocated. Further, hi = ih and /3j = P,i = 
1, 2, 3, 4 where j < P < |. Any profile in which two mobiles 
associate with one BS, and the remaining three with another 
is a NE. On the other hand, (ai = a.2 = 1, 03 = an = as = 2) 
and (ai = <22 = 2, 03 = 04 = 05 = 1) are the only socially 
optimal NEs. 

3) Collocated Mobiles, Symmetrically Placed BSs, and Sin- 
gle Class Traffic: In the special case when all the mobiles 
are collocated and all the BSs are symmetrically placed with 
respect to the collocated mobiles, we have the following result. 

Proposition 3.6: All the NEs in the game (A4,Af, (cj,i G 
with h^ = h for all i G A4,j G Af, are system optimal. 
Proof: The mobiles as well as BSs are indistinguishable 
in this game. At a NE, let rrij be the number of mobiles 
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associated with BS j. We first prove that at any NE, the vector 
of mobiles' costs is unique up to permutations. To prove this, 
it suffices to prove that the vector m = (rrij , j £ Af) for a NE 
is unique up to permutations. As m yields a NE, the following 
must hold for all j, k G Af: 



P 



< 



hi — nijf3 hi — rrikfi — f3 

< m k + 1. 



3 ■> 

or m , 



(4) 



Define n — and I = M — nN. From (O we see that m 
given by rrij = n + 1, j — 1, . . . , I, rrij = n,j = I + 1, . . . , N 
characterizes one of the NEs; other NEs are permutations of 
this vector, and m is unique up to permutations. We now show 
that m is a system optimal congestion vector, and the system 
optimality of all other NEs follows. To do this observe that 



c(a) = t £ r 



j3 a 1 rrij (5 



is a Schur-convex function in (mi, . . . , mjv) because is 
a convex function. This implies that the minimum value is 
attained at a vector which is as close to uniform as possible, 
i.e., a vector that is majorized by any other vector (Marshall 
& Olkin (23)0 

All such vectors are permutations of m (Al- 
ternatively, if there exist BSs j and k such that rrij > m k + 2, 
moving a mobile from BS j to BS k results in a strictly lower 
cost). This concludes the proof. ■ 

IV. Continuum of Mobiles 

In this section, we consider a nonatomic version of the 
system in Section IIII-AI Such a model is of interest for 
two reasons. First, for many of the fixed QoS traffic classes 
(e.g., voice), the target SINR requirements in CDMA cellular 
systems are very small. In a typical IS 95 CDMA system with 
system bandwidth 1.25 MHz, chip rate 1.2288 Mcps, data rate 
9.6 Kbps, and target || = 6 dB, the target SINR turns out to 
be -15 dB, i.e., (Kumar et al. 11251 Chapter 5]). If we assume 
that at any time the number of mobiles associated with a BS 
is large, it is reasonable to say that an incoming mobile or 
an outgoing mobile has a negligible effect on the congestion. 
Secondly, we have seen that our proposed algorithm may end 
up with inefficient associations. There is extensive work on 
toll mechanisms that induce system optimality in networks 
with a continuum of users. The analysis of toll-mechanisms (or 
pricing) on a multichannel multicell network with a continuum 
of mobiles can be expected to shed light on the existence and 
properties of pricing mechanisms for networks with discrete 
mobiles. 

A. System Model 

Let M — ufL-^Mi be an infinite set of £ = {1, . . . , L} 
classes of nonatomic mobiles. By nonatomic mobiles, we 
mean that the effect of a single mobile at a BS is infinitesimal. 
The population of class I mobiles has "mass" Mi. All the 

7 The condition hij = hi for all j 6 J\f is used to deduce that NE profiles 
ny i 

is used to deduce Schur-convexity of C(a). 



mobiles in a class are collocated and require equal minimum 
SINR. In particular, all such mobiles have same power gains 
to any of the BSs (gains from a mobile to different BSs can be 
different). Assume Af to be the finite set of BSs. As before, a 
denotes the common standard deviation of receiver noise at all 
BSs. Let 7; be the common minimum required SINR density 
for class I mobiles, and hij be the power gain between a class 
I mobile and BS j. An association profile a is a measurable 
function a : A4 — > Af. Any association a leads to a congestion 
profile m(a) = (mij(a),l € C,j € Af), mij(a) being the 
mass of class I mobiles associated with BS j. Let A4* denote 
the set of all such congestion profiles. 

Under an association profile a and a power density allo- 
cation p : M ->• K + , the SINR density for x € Mi, I € C 
is 

hl a ( x )p(x) 



Ezti S M , z)hi a[z) p{z)dz + a 2 ' 



where 



l a (x,z) 



1 . if a(x) = a(z) 
0. otherwise 



Our definition of a "class" makes all the mobiles in a class 
alike, and so, congestion profiles are sufficient to characterize 
the system. In the sequel, we just use mij for mij(a) for 
convenience. The dependence on a is understood. 

Consider again the subproblem of power control with a fixed 
congestion profile m. The following result is analogous to 
Proposition 13.11 and is shown in Appendix lAl 

Proposition 4.1: 1) The power control subproblem of BS 

j is feasible iff J2iec m ij^li < 1- 
2) If the power control subproblem of BS j is feasible, there 
exists a unique Pareto efficient power density p given by 

p(x) 
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Vx G AAi such that a(x) = j,l e C, where a is the 

underlying association profile. 
An evolutionary dynamics can be proposed to address the 
combined association and power control problem. To this end, 
we define functions cy : A4* —> R+, where cy(m) denotes 
the minimum power density for class I mobiles associated with 
BS j under congestion profile m, as 

c ,x = 7^! 

For notational convenience, define 



9lj 



and 



We then have 



c(z) 



Clj( 



hij ' 

L 

= y^7;mtj,Vj E Af 
i=i 

= f T i iJ if*<l 
|_ oo, if z > 1. 

= 9ijc(rrij). 



(5) 



Again we assume that the system is feasible, i.e., there exists 
a feasible assignment, as done in Section IIII-BI This boils 
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down to assuming J^iec ll-Mi < N in the case of nonatomic 
mobiles. Now, structures of the cost functions allow us to 
restrict attention to the region where rrij < l,Vj G Af; if 
rrij > 1 for a j G A/", all the mobiles associated with j incur 
infinite cost. 



B. A Congestion Game Formulation 

We model the problem as a nonatomic congestion game. 
The continuum of mobiles constitute the population, and Af 
denotes the common action set for players of all the classes. 
Class I players are accompanied by cost functions cy(m),j E 
Af. In the following, we refer to it as the game (A4 , Af, (cy , I E 

£,j e A/")). 

Proposition 4.2: The nonatomic game (A4,Af,(cij,l E 
£,j E Af)) is a potential game. Furthermore, it admits at 
least one NE, and the set of NEs coincides with the set of 
minimizers of the potential function. 

Proof: In the region {m : rrij < 1, Vj E Af}, the function 
V : M* -> E U {+00} defined as 



V(m) := ^ [X^TiWy loggy + / 
jeM \iec Jo 



log c(x)dx J , (6) 



is a C 1 function with 
dV(m 



dm 



li log Sy + 7; log c{m 3 ) = 7, log cy (m) 



for all £ € £, j € A/", m E A^*. Thus the nonatomic game 
(M , A", (7; log cy , / € £,j G A/")) is a potential game with 
F(m) as a potential function (see Definition 12.51 . Note that 
the strategic form game (A4,N , (cy, I E £,j E Af)) is better 
response equivalent to (A4,Af, (7; logcy, / € £,j G A/ - )). 
Thus the former is also a potential game with the same 
potential function V(m). 

Now consider the following optimization problem 



Minimize ^( m ) 

subject to >J my = Mi, I e £ 

mij > 0, I E £,j E Af. 
All the conditions are self-explanatory. Observe that 



(7a) 
(7b) 



<9 2 l/(m) 



7 l 7zc(™ i ) 




if fc = j, 
otherwise. 



drriikmij 

Thus, with an appropriate ordering of the components of m, 
the hessian of V^(m) is given by 

"c(mi)D ... 
c(m 2 )D ... 



V 2 F(m) = 











where 



L> := 
and T := 





rr T 

[71, • 



c(tojv)D 



(8a) 
(8b) 



Clearly, D, and hence V 2 V^(m) is a positive semi-definite 
matrix. Thus, V(m) is a convex function of m. Since we 



are minimizing a convex objective function subject to linear 
constraints, there exists at least one minimizer, and all minima 
are global minima. Also, Kuhn-Tucker first order conditions 
are necessary and sufficient 11281 Section 5.5.3]. Combining 
this with the fact that NEs are the profiles which satisfy the 
Kuhn-Tucker first order conditions for a minimizer of the 
potential function (see Section lll-Bb . we see that the set of NEs 
coincides with the set of minimizers of the potential function. 

■ 

Remark 4.1 : The assertion in the above proposition does 
not hold for general population games. While all local max- 
imizers of potential function are equilibria, not all equilibria 
maximize potential (even locally) in general lfT7l Section 3]. 
This is unlike finite player potential games where only equi- 
libria are the local maximizers of potential functions. 

Furthermore, NEs have the following property (see l29l 
Proposition 3.3]). 

Proposition 4.3: The congestion at a BS is constant across 
all the NEs of the game (M,Af,(ctj,l E £,j E Af)). 
Consequently, the cost density for a class is also constant 
across all the NEs. 

Remark 4.2: At NEs, the congestions (at BSs) by class, 
my, are not unique because the objective function <j6j is not 
strictly convex with respect to this set of variables. 

C. System Optimality 

Analogous to the definition in Section IIII-DI we define the 
system performance measure 



C(m) 



L 

EE 

3 eM 1=1 



mij9ijc(mj). 



(9) 



A congestion profile m* E A4* is said to be system optimal 
if it minimizes C(m) over all possible profiles m E A4* . 

In contrast with the discrete mobiles case where equilibria 
need not be Pareto efficient (see Example 13. U . we have the 
following result for the nonatomic case. 

Proposition 4.4: All NEs of the nonatomic game 
(M,Af, (cy, I G £, j G AT)) are Pareto efficient. 

Proof: Let m be a NE congestion profile. Under a NE, 
the cost densities for the mobiles of the same class are equal, 
irrespective of their associations (see Remark l2~2l i. Thus, it is 
sufficient to prove that there does not exist another congestion 
profile m' such that for every class and for all BSs j, fc, 
with mij > 0, m' lk > 0, 



Cifc(m') < cy(m), 



(10) 



and strict inequality holds for some such l,j and k. Assume 
that such an m' exists. Then, 

9ikc(m' k ) < gijc(mj) < gi k c{m k ) 

where the last inequality follows because m is a NE and mij > 
0. This yields m' k < m k . This further implies that there is a 
BS s such that m' s > m s , and a class t such that m' ts > m ts . 
By the strictly increasing property of c, we have 

gt s c(m' s ) > g ts c(m s ) > g tr c(m r ) 
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for a BS r such that m tr > 0. Such a BS of course exists 
and the latter inequality follows because m is a NE. The two 
inequalities imply Cf S (m') > c tr (m), and so the tuple t,r,s 
violates ( fTOb . Thus the assumption that m' Pareto dominates 
m is incorrect. This completes the proof. ■ 
We show that the NEs are system optimal if all the mobiles 
are collocated, and all the BSs are symmetrically placed 
around them. 

Proposition 4.5: All NEs in the nonatomic game 
(M,J\f,(cij,l G C,j G TV)), with hij = h for all 
/ G C,j G Af, are system optimal. 

Proof: In the case of collocated base stations 

a 2 L 
C( m ) = -r^E w 3 c W 

j€Af 1=1 
= m A m i) 

For the reason described earlier, we restrict attention to the 
region where rrij < 1 , Vj G Af. In this region, 



-m ]C (m 3 ) 



1 



dm j (1 — mj 



and so mjc(m,j) is a convex function of rrij. Thus C(m) is a 
Schur-convex function of (m-,-, 1 < j < N), and is minimized 
at any m* with 



TTLa 



lec 



for all j G TV. When fry = ft for all / G C, j G Af, any 
congestion profile with equal congestion at all the BSs is a 
NE. Thus, the system optimal profile m* is a NE. Since all 
the NEs incur equal cost (see Proposition 14.3b . all NEs are 
system optimal. ■ 

However, NEs need not be system optimal if BSs are not 
collocated, or mobile are not collocated. We illustrate these 
facts through the following examples. 

Example 4.1: Consider an infinite set AA of nonatomic 
mobiles belonging to two classes; class 1 and class 2 mobiles 
have masses M\ and 3Mi respectively. Assume common 
minimum SINR density requirement 7, and let 3M17 < 1. 
Let there be two collocated BSs. Let the power gain between 
a class I mobile and a BS be hi, hi < -y-. A congestion profile 
is a NE if and only if it assigns equal load to both the BSs. 
Thus, the total cost incurred at NE 



C* 



+ 



hi(l-2jMi) /i 2 (l-2 7 Mi) 

Next, consider a profile in which class 1 mobiles associate 
with BS 1 and class 2 mobiles associate with BS 2. The cost 
incurred now is 

1 M 1 a 2 3 7 M 1( t 2 
fci(l-7Afi) fca(l-37Mi)' 
It can be easily checked that C < C* if 

1 h 2 /'i - h x 



Mi < 



7 



hi 



Example 4.2: Consider an infinite set AA of nonatomic 
mobiles all belonging to same class; M := \M\. Assume com- 
mon minimum SINR density requirement 7, and let A/7 < 1. 
Let there be two BSs with hj the gain to BS j, j = 1,2. An 
NE congestion profile (a*M, (1 — a*)M) is given as 
1) if < (1 - Mi), a* = 0, 



2) if | < (1 - Mj), a* 

3) otherwise, a* satisfies 

7<t 



1, 



7<7 



hi(l - a*7M) 
1 - a*jM 
1 - (1 - a*) 7 M 



Mi-(i 

h2 

hi 



*) 7 M) 



(11) 



On the other hand, a congestion profile (a°M, (1 — a°M) 
will be system optimal if and only if a° solves the following 
optimization problem: 



ajMa 2 



Minimize 

hi(l - cryM) 

subject to < a < 1. 



(1 - a^AIa 2 
h 2 {l - {1 - a)jM) 



(12) 



This is a convex optimization problem, and it is straightfor- 
ward to show that 

- A/7), a° = 0, 



2) if < (1 - Mi), a° = 1, 

3) otherwise, a° satisfies 

1 - a°jM 




1 - (1 - a°)jM 

Hence, if min{^, > 1 - A/7, then minjj^, \j j^} > 
1 — A/7, and a* and a° must satisfy ( fTTT > and ( TT3b respectively. 
In such a case, the NE will be system optimal if and only if 
hi = h 2 . 

Remark 4.3: Sandholm ifTTl shows that if the cost function 
for each mobile is a homogeneous function of a certain degree, 
then all NEs are system optimal. Note that in Example 14.21 
NEs are not system optimal unless hi = h%. We remark that 
the system optimality for the latter case does not follow from 
Sandholm |17| because the cost functions are not homoge- 
neous functions. 

V. Pricing for System Optimality 
A. Continuum of Mobiles 

In this section, we show that there is a toll mechanism that 
can induce system optimal associations and power allocations 
in a cellular network with multiple classes of mobiles. We also 
show that the mechanism can be employed in a distributed 
fashion. 

Define 



c'(z) 



J dz 

00 



if z < 1 
if z > 1 



Consider a congestion profile m = (mij,l G C,j G Af). We 
propose that a class I mobile joining BS j be levied a toll 

L 

(14) 



^■(m) = -1^, „,,</, ,<'{„>,). 
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Now, define cy(-) = cy(-) + %(•)? Vi G £,j G TV, and 
consider the nonatomic game (A4, Af, (cij, I G £,j G A/ - )). 
Players may incur different power costs (cy(-)) in different 
NEs of this game. Therefore, one has to distinguish between 
the following two cases (see Fotakis & Spirakis lf30l ). 

1) A toll mechanism is said to weakly enforce system 
optimality if some NE of the game with tolls is an optimal 
profile. 

2) It is said to strongly enforce system optimality if all the 
NEs of the game with tolls are optimal profiles. 

We show that tolls %(•) weakly enforce system optimality in 
all cases and strongly enforce it in a special setting. 

Proposition 5.1: The nonatomic game (A4,Af,(cij,l € 
C,j G Af)) is a potential game. Furthermore, a congestion 
profile m is system optimal only if it is a NE of this game. 

Proof: Recall the system performance measure C(m) 
defined in (O. Observe that 

dC(m) _ 

for all I G C,j G Af, and m in {m : rrij < G Af}. Thus 
the nonatomic game (A4,Af, (cy, I G C, j G Af)) is a potential 
game with C(m) as a potential function (see Definition 12.51 ). 

Next, observe that system optimal associations and powers 
are solutions of the following nonlinear optimization problem: 

Minimize C( m ) 

subject to Conditions ( TTab - ( 1751 ). 

Since all the constraints are linear, any optimizing congestion 
profile of C(m) necessarily satisfies the Karush-Kuhn-Tucker 
first order conditions l28l Chapter 5]. But, any congestion 
profile satisfying these conditions is a NE of the game 
(A4,Af, (cy, I G C,j G Af)). Thus, any system optimal 
congestion profile is also a NE of this game. ■ 

If all the mobiles are collocated, the proposed tolls strongly 
enforce system optimality. 

Proposition 5.2: All NEs in the nonatomic game 
(M,Af,(cij,l G C,j G Af)), with hij = hj for all 
I G C,j G Af, are system optimal. 

Proof: It suffices to show that C(m) is a convex function 
if ft,y = hj for all i G C, j G TV. Then, any congestion profile 
satisfying Karush-Kuhn-Tucker conditions (i.e., any NE) is 
system optimal. 

If fly = hj for all I G C, j G /V, 

<7(m) = £ 

Using the observation c(a;) + xc'(x) = c'(x), it is easy to see 
that 

d 2 C(m) = j 7 ' 7 '\ ( " j)ct2 if fc = j, 
dm ik mij I otherwise, 



and 



V 2 C*(m) 



e"(mi)<T 2 £) n 



o 






c"(m N )a 2 



D 



with £> given by (f8ab-(t8bT>- It is now obvious that V 2 C(m) 
is a positive semi-definite matrix, and so C(m) is a convex 
function of m. ■ 
However, the tolls ty(m) may fail to strongly enforce a 
system optimal congestion profile even if all the BSs collo- 
cated, the mobiles require a constant SINR density 7, but they 
are not collocated. To see this, consider the congestion profile 
m* with 

m y = ^Vle£,j eAf. 

It can be easily checked that, for all I G C, 



cy(i 



7cr 2 /7M 



E 



7 2 a 2 Mi 



,h_M_ 



N \ N 



which is independent of j G A/". Thus m* is a NE of the game 
(A4 , Af, (cij , / G C,j £ Af)). But m* may not be system 
optimal (see Example 14. U . 

Remark 5.1: 1) cy = Cy + iy can be interpreted as the 
marginal cost due to additional association of class / mobiles to 
BS j. The term cy is the power density incurred by these new 
mobiles, and iy is the increase in power consumption densities 
of the mobiles already associated with BS j, integrated over all 
such mobiles. Economists call them "private cost" and "social 
cost", respectively. Selfish mobiles do not care for the social 
cost, while the social optimality criterion accounts for this 
marginal externality BP . 

2) The cost functions for various classes have a certain 
structure in the settings of interest to us. Mobile classes that 
consider a BS pay tolls proportional to their required SINR 
densities. In particular, tolls are uniform across all the mobile 
classes that have equal SINR requirements. This is special to 
our setting; usually one does not see uniform tolls in the case 
of multiclass networks (see Dafermos |20|, Smith [21]). 

This toll mechanism can be implemented in a distributed 
fashion. All the BSs broadcast the tolls (normalized by SINR 
densities) along with their aggregate congestions as before@ 
All mobiles need to know their scaled gains to each BS 
j G Af. A mobile then makes a choice taking both power 
density and toll into account. 



B. Discrete Mobiles 

Pricing mechanisms for networks with discrete mobiles 
are relatively difficult to design and analyze (Fotakis & Spi- 
rakis l30l ). Again, we propose a toll mechanism that weakly 
enforces system optimality in all cases and strongly enforces 
it in a special setting. The mechanism is motivated by the toll 
mechanism for the nonatomic case (Theorem 15. It . 

Consider the network model of Section IIII-AI and an asso- 
ciation profile a'. Let mobile i evaluate BS j for association. 
Define a = (j, a'_ ? ). Analogous to the nonatomic case, define 



8 Normalized tolls are uniform across all mobile classes that consider 

11 

a BS. A mobile can recover the exact toll from the normalized value. 
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"private" and "social" costs as 



Cj(a) 



Pi 



h ij I 1 - EfcsA-l^a) ftkY 



ft, 



andti(a) = V _ — 

ltM%\{i} hlj ^m 3 (*) ftk\ + 



ft, 



(15) 



I 1 EfceM J (a)\{i} ftk] + 

respectively]^ Clearly, Ci(a) is the required power of mobile 
i if it joins BS j, while U(a) is the aggregate increase in 
power consumption of all other mobiles associated with BS j. 



M 



given 



We propose a toll mechanism with tolls ti : Af 
by (15[ . This yields a new game (M, Af, (gj, i G M)) with 
cost functions for an association profile a given by 

c j; (a) = Ci(a)+ti(a) 

= y £ 



hi* [1-E 



£ £ 



(16) 



Proposition 5.3: The finite strategic form game 
(M,Af, (ci,i G A^)) is an ordinal potential game and 
thus admits FBRP. 

Proof: For the game (M,Af, (c-i,i G M)), the function 
F : A/" |A11 -> R given by 

A 



is an ordinal potential function, as can be straightforwardly 
checked. Thus (M,Af,(ci,i G M)) is an ordinal potential 
game. Since it is also a finite game, the FBRP property holds. 

■ 

Note that the potential function V(a) equals the system 
performance measure C(a) defined in Section IIII-DI Hence 
an association profile a° that optimizes system performance 
is also a (global) minimizer of V(a), and therefore a NE of 
the potential game with tolls. So, we see that tolls tj(a) weakly 
enforce a system optimal association profile. In general, tolls 
do not strongly enforce a system optimal association profile. 
For instance reconsider Example 13.11 The association profile 
(di = 2,a2 = 1) is inefficient, but an NE for the game 
{M,Af, faieM)). 

In the following we consider special cases, and investigate 
the effect of the proposed tolls. 

1) Collocated Mobiles with Single Class Traffic: Let us 
consider the special case when all the mobiles are collocated 
and have identical minimum SINR requirements. In other 
words, hij = hj and fti = ft for all i <G M,j <E Af. The 

9 In 115t . when both terms within parentheses are oo, the expression is taken 
to be oo; we may think of driving /3 to the true values from below, and the 
first term always dominates the second. Same remark holds for other such 
expressions also. 



potential function for this special case can be written as 

|JWi(a)|/J 



hj [1 - \MM)\ft] 



Define 9j = g,/(m) = and m 3 (a) = |^-(a)| for 

all j G Af. Then m(a) = (rrij(a),j G J\f) denotes the con- 
gestion profile under a. Since mobiles are indistinguishable, 
any two association profiles that lead to identical congestion 
profiles are essentially indifferent from the point of view of 
analysis. Thus we talk solely in terms of congestion profiles. 
Abusing notation (the argument of V(-) was earlier defined to 
be the association profile a), we write 



V(m) 



Since (M, Af, (cj, i G M)) is a finite potential game, an 
association profile m* will be a NE if and only if 

9jf(m*) + g k f(m%) < 9j f(m* - 1) + g k f(m% + 1) (17) 

for all k j, j, k G Af. The following proposition shows 
that tolls tj(a) strongly enforce a system optimal association 
profile in case of collocated mobiles with single class traffic. 

Proposition 5.4: All the NEs in the game [M , Af, i G 
M)), with hij = hj and ft t = ft for all i G M,j G TV, 
are system optimal. In other words, the tolls strongly enforce 
system optimality. 

Proof: Let m° be a system optimal congestion profile, 
and m* any other profile such that V(m*) > V(m°). Partition 
the set Af as Af = Afo U Af+ U Af- such that 



j G Afo 
j G Af+ 

j G Af- 



m* > m° + 1 



to* < m°j 



1 



Start with the congestion profile m*, and move mobiles 
from BSs Af+ to BSs Af- one mobile at a time, so that we 
end up with the congestion profile m°. In this process we get 
a succession of congestion profiles, each of which satisfies 

rrij = m* V j G Afo 
rrij < m* V j G Af+ 
rrij > m* V j G Af- 

There must exist a pair of successive congestion profiles 
m' and m" such that V(m') > V(m"), with m" possibly 
the ultimate congestion profile m°. Let m" be obtained from 
m' by the transfer of a mobile from BS j G Af+ to a BS 
k G Af- . We then have 

9jf(m'j) + g k f(m' k ) > gjfim'j - 1) + g k f(m' k + 1) 

which is same as 

9j(f(m'j) - f(m'j - 1)) > g k (f(m' k + 1) - /(m' fe )). (18) 

Recall that / is a convex function and m'j < m*,m' k > m* k . 
Using these in (fT8l , we get 

ft(/K) - f(m* - 1)) > g k (f(m* k + 1) - f(m%)), 
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i.e., 

9jf{m*) + 9kf(m* k ) > g 3 f{m* - 1) + g k f{m* k + 1) 

which implies that m* is not a NE (see (T7\). This completes 
the proof. ■ 

2) Collocated Mobiles and Symmetrically Placed BSs: 
Now we consider another special case when all the mobiles 
are collocated and all the BSs are symmetrically placed with 
respect to the collocated mobiles. In this case hij = h for all 
i S M,j € Af. We have the following result. 

Proposition 5.5: With hij — h for all i £ M,j £ Af, the 
NEs in the game (A4,Af, (ci,i £ .M)) coincide with those in 
{M,Af, (ci,i€M)). 

Proof: With hij = h for all i G M , j £ TV, 

Pi 



Cj(a) = 



h [l-E ;e M ai (a)A] + - 



Thus an association profile a is a NE of the game 
(M,Af, (ci,i € Al)) if and only if 



E 

for all j £ Af \ {a.;} and i £ M. Next, 
c 2 / J2ieM a .{a) Pi 



(19) 



c,(a) 



2~2l£M at (n) Pi ~ Pi 
I 1 - Sie7W Q! (a)\{i} A]" 



h [1 - E; e ^ a ,(a) AH 1 - EieA^MU*} A] + ' 

It can be easily checked that NEs of the game (M,Af, (c«, i £ 
M)) are also characterized by (fl~9b . Hence we get the claim. 

■ 

Thus tolls may not strongly enforce a system optimal associ- 
ation profile in this case (see Example 13.31 ). 

3) Collocated BSs with Single Class Traffic: Even in this 
special case tolls tj(a) may fail to strongly enforce a sys- 
tem optimal association profile. For an illustration reconsider 
Example 13.41 The association profile (<zi = 03 = 1,0,2 = 
124 = 05 = 2) is not system optimal, but an NE for the game 
(M,Af, (5i,i£M)). 

Remark 5.2: 1) While tolls at a BS are equal for all 
the mobiles not associated with it and having equal SINR 
requirements, they are mobile dependent for all associated 
ones (see (TTBTl). This is unlike in nonatomic case where we 
saw uniform tolls at a BS for all the mobiles with equal SINR 
requirements. 

2) The modified algorithm (the one accounting for tolls) can 
be implemented in distributed fashion. All the BSs broadcast 
quantities t° (a) given by 

5W- E T. * 



Ze7Wj(a) 



hij [1-E 



k£M j (a) 



Pk] 



along with their aggregate congestions Efeejw -(a) Pk- the 
mobiles need to know the scaled gains -4f of their own 



channels to all the BSs j £ Af. Mobiles use these broadcast 
information to calculate their powers and tolls, and choose a 
BS taking both into account. 

Discussion: The proposed pricing technique can be used 
to induce a system optimal routing in atomic weighted network 
congestion games with arbitrary nondecreasing edge latency 
functions |f32|F"l In this setting, the joint BS association and 
power control problems can be viewed as network congestion 
games over two-terminal parallel-edge networks: the edges are 
identified with BSs, and latencies are identified with minimum 
power requirements. It turns out that the proposed tolls weakly 
enforce a system optimal routing profile in general network 
congestion games. They strongly enforce a system optimal 
routing profile if 

1) the network is two-terminal series parallel, 

2) the users are unweighted (i.e, have identical weights), and 

3) the latency functions are standardFl 

VI. Conclusion 

We studied the combined association and power control 
problem in multichannel multicell cellular networks. We stud- 
ied the cases of discrete mobiles and a continuum of mobiles. 
We proposed several distributed mechanisms motivated by the 
techniques of game theory. We studied the inefficiency of the 
distributed algorithm in the case of a continuum of mobiles. It 
is an open question whether such inefficiency can be quantified 
in the case of discrete mobiles. To mitigate the inefficiency, 
we proposed toll mechanisms in both the settings. 

Throughout we assumed only one BS per channel. It would 
be of interest to extend the work to multiple BSs per channel. 

Appendix A 
Nonatomic Power Control 

Assume M to be an infinite set of mobiles, and a single BS. 
The required SINRs of mobiles are given by the function 7 : 
M — > Let h : M —> K++ be the gains of mobiles and 

let p : M — >• K ++ be the power allocation. The functions p(x) 
and 7 (a;) are interpreted as power density and target SINR 
density, respectively, per unit mass. The feasibility condition 
for p can be written as 

p(x)h{x) 



I M p(y) h (y) d y + °' 

or equivalently, 



> j(x), Vx 6 M, 



p{x)h(x) > j(x) / p(y)h(y)dy + r y(x)a , Vx £ M. 
J M 

Integrating the inequalities over the set of all mobiles, we get 
the following necessary condition for feasibility. 



p{x)h{x)dx > / j(x)dx / p(y)h(y)dy 
Jm Jm 

+ a 2 [ j(x)dx (20) 
Jm 



> 



7(x)dx / p(y)h(y)dy 
m Jm 



10 Here, the system cost is weighted sum of the latencies of all the users. 
"A latency function c( ) is called standard if mc(m) is convex 1311 . e.g., 



c(m) = -r^— 

v ' 1 — 7] 
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where the strict inequality arises because a 2 J M r y(x)dx > 0. A. Two BSs 



Canceling J M p(y)h(y)dy from both sides, we get the neces- 
sary condition 



j(x)dx < 1. 



M 



Assuming ( l2Tb holds, 



p{x) 



"/(x)a 2 



h(x)(l - J M j(x)dx) 



(21) 



(22) 



is a feasible power allocation, and so (TZTT i is necessary and 
sufficient for feasibility of the power control problem. 

The power allocation given by d22l is Pareto efficient. Let 
us see why. Suppose q : A4 — > K++ is another feasible power 
allocation such that q(x) < p(x) and strict inequality holds 
for a set of mobiles having a positive measure, and so 

q{x)h(x)dx < / p(x)h(x)dx, 
M JM 

a ' 2 Sm li x ) dx 

1 - Sm li x ) dx 

where the equality uses (l22l . Thus, q violates the necessary 
condition 

I \hi \j >» g2 Im l( x ) dx 
q(x)h(x)dx > p — r ~r- r , 

m 1 - J M Jl x )dx 

a rearrangement of d20l i, and hence cannot be feasible. 

In fact p can be shown to be the unique Pareto efficient and 
hence the system optimal power allocation^! Indeed if q is 
another feasible power vector which is also Pareto efficient, 
the pointwise minimizer of p and q is also feasible and 
Pareto dominates p, thus contradicting the fact that p is Pareto 
efficient. 



Appendix B 
Price of Anarchy: Continuum of Mobiles 

Recall that a NE is not necessarily a system optimal con- 
gestion profile (see Example 14.21 ). Coordination ratio IT331 or 
Price of Anarchy l34l characterizes the inefficiency caused by 
the selfish behavior of players; it is the ratio of the cost of the 
worst NE and the optimal cost. We observed in Proposition l4.3l 
that, in the nonatomic case, mobiles incur the same cost at all 
the NEs. We can then define price of anarchy as follows. 

Definition B.l: Let m be a NE, and m° be a system optimal 
congestion profile. Then the price of anarchy is 



PoA 



C(m) 
C(m°)' 



We restrict our analysis to a single class population. We 
assume that all the mobiles have identical minimum required 
SINR density 7 and identical power gain hj to BS j, j E Af. 



The system optimality criterion is to minimize the sum of power con- 
sumptions over all the mobiles. 



First we consider a case with 2 BSs as in Example 14.21 
Let hi > /7 2 Also, let (a*M, (1 - a*)M) and (a°M, (1 - 
a°)M) be the congestion profiles under a NE and a system 
optimal association, respectively. Recall from Example l4.2l that 

1) if 7M < 1 - Jfe, then a* = a° = 1 

2 ) if 1 - v/if < 7 A/ ^ 1 - If > then a * = !> and from CL! 
\fh\ - \fh~2 + ^M\fh~2, 



a 



3) if 7 A/ > 1 



(Vhi + VfohM 
then from (QT) 

hi — hi + 7M/12 



a 



(/n + /i 2 ) 7 M 
and a° is as above. 
C(m) and C(m°) are obtained via substituting a = a* 
and a = a", respectively, in the objective function (f72l) . 
Straightforward calculations give that 

if M < 

A(2 — r M)fM 



PoA(Af) 



where A 



1 



_ } (l- 7 M)(2v / A-(l- 7 M)(l + A)) 
7A/(1+A) 



if < M < ^ A 



L 2v / A-(l-7M)(l+A) 



7 — — 7 

if M > ±=± 

— 7 



< 1. Further calculations also yield that 



PoA(M) is continuous at M 



= (l-A) 
7 ' 



and 



dPoA(M) J > if M < 



dAf 



< if M > ±=± 

— 7 



Thus, the price of anarchy is maximized when M = 
Moreover, the maximum price of anarchy is 

1 - A 2 



l-A 



2\/A - A(l + A) 

Viewing this now as a function of A G (0,1], we see that the 
maximum price of anarchy decreases with A. We also observe 
that PoA — > 00 as A — > 0, i.e., arbitrarily high PoAs can be 
realized in 2 BS networks. 

B. N BSs 

Again, without any loss of generality, we assume that 
hi > h,2 > ■ ■ ■ > hnr- We also assume that the population's 
mass is Aj when it spills over BS j under NE. Clearly, 
A2 < A3 < • • ■ < A at. In the case of 2 BSs we proved 
that price of anarchy is maximized when the population spills 
over BS 2 under NE. In the case of N > 2 BSs also, 
simulations suggest that the price of anarchy is maximized 
at one of the spill over points {Aj, j = 2, . . . , N}. We have 
however not been able to prove this observation. We illustrate 
this observation in Figure [T] which shows the price of anarchy 
as a function of mass of the population in a network with 
5 BSs. For this example 7 = 0.01, and the power gains are 
chosen independently and uniformly from [0, 1]. 

13 If hi = h,2 equal fraction of population join each of the BSs under the 
NE and the system optimal association, and the price of anarchy is 1. 
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Fig. 1 . PoA vs M : the price of anarchy is maximized at one of the masses 
{A 2 ,A 3 , A 4 , A 5 }. 



In the following, we give expression for the price of anarchy 
when the total population has mass Aj,j > 2. Let mj, k < j, 
be the mass of the population associated with BS k under NE. 
From the definition of NE (see Definition \2.4j . 



7CT 



for all k < j. 



Thus, 



where ej 



hk(l-jm k ) hj 

J2 k<j for j > 1. The cost at NE 

n< \ ^ 2 f J ~ 1 

C(m)=^m,— = a f — - 

k<j 3 v 3 

Next, consider the system optimal congestion profile m* . It 
can be easily checked that m k is positive only if k < j 1 for a 
j < j' < N. Moreover, it must be the case that 

\/hk(l - ym%) = K for all k < j', 
and Jhji < K, 

for some constant K > 0. Recall that J2k<j' m *k = Ar 
Straightforward calculations yield that 

/ -j + hjej-i 



K = 



-l-i 



where e* := X)fc<j f° r i — ^ Moreover, j' is the least 
index satisfying 



-l-l 



The optimal cost, after substituting the values of m* k , turns out 
to be 

C(m*) = J2 ml 



7(7 



fc<j' 

A 



'■hk{\-m* k i) 



e f-i 



ej'-i 



The price of anarchy PoA(Aj) is the ratio §^0k ■ 

Finally, we show that the price of anarchy decreases with 
mass for M > Ajv- It can be easily checked that, for M > 

Ajv, 



PoA(M) 



e N Mj~ (e N N - et) 



1 



e N N 



' N 



e N Mj - {e N N - e* N ) 



from which the claim follows. Thus, to obtain a bound on the 
price of anarchy, we only focus on M < Ajy. For M < Ajv, 
the load on BS j 

rrij < - 1 - — 
7 V hj 

under NE. We use this observation in the next section. 

C. A Bound on the Price of Anarchy 

Now, we derive a sharp bound on the price of anarchy 
for single class networks with arbitrary number of BSs, 
and gains hj G [h 
Roughgarden [31] Chapter 3] 

In the BS association game, a generic cost function is of 
the form 



„, /i max ] for all the BSs. We follow 



c h (m) :-- 



1 



h 1-jM' 



and 



C := {Ch{-) : h € [/l m in,/lma 



]} 



is the class of all feasible cost functions. Observe that the 
functions Ch(-) and the class C both are standard^ We define 

d(mch(m) 



c h (m) := 



dm 



We also assume that the load on a BS with gain hj does not 
exceed 

7 V h 

under NE. Thus, we redefine anarchy value for a cost function 
c h (-) a£! 

a(ch) := sup [A/x + (1 - A)]" 1 

m<6 h 

where A <G (0, 1) satisfies c/,(Am) 



c/j (m) and fi 



Ch < 1, Both A and u are functions of to; we do not 
show this dependence explicitly. Straightforward calculations 
yield that 



A 
/' 

and a(ch) 



1 — \/l — TO7 
TO7 



my, 




14 A cost function c(-) is called standard if mc(m) is convex. A class C i s 
standard if it contains a nonzero function and if each c(-) £ C is standard 1311 . 

15 The original definition 1311 Definition 3.3.2] considers supremum over 
m S (0, oo). 
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The anarchy value for class C is (see 13T1 Definition 3.3.3]) 

a(C) = sup a{c h ) = \ ( 1 + t/^ 1 ] 

c h GC ^ \ V thorn I 

It can be easily checked that OTl Theorem 3.3.8] remains 
valid with our new definition of anarchy value. Thus, price of 
anarchy is bounded by a (C ) . For any < e < a (C ) — 1 , a price 
of anarchy > a(C) — e is realized in a network in which (i) 
there is one BS with gain /i max , (ii) there are several BSs with 
gain h m - m (minimum number depending on e), and (Hi) the 
population has mass 6'/ tmax (see the proof of [[3T] Lemma 3.4.3] 
for details). 
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